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HIGHLIGHTS 


►  Redox  flow  batteries  coordinate  with  IPFC  on  Load  Frequency  Control  investigated. 

►  BFO  algorithm  is  used  to  optimizing  the  controller  design  parameters. 

►  Control  Performance  Standard  based  controller  reduction  wear  and  tear  of  equipment. 

►  Simulation  result  gives  coordinate  control  improves  dynamic  quality  of  Power  system. 
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This  paper  proposes  a  sophisticate  application  of  redox  flow  batteries  (RFB)  coordinate  with  Interline 
Power  Flow  Controller  (IPFC)  for  the  improvement  of  Load  Frequency  Control  (LFC)  of  a  multiple  units 
two-  area  power  system.  The  Interline  Power  Flow  Controller  is  to  stabilize  the  frequency  oscillations  of 
the  inter-area  mode  in  the  interconnected  power  system  by  the  dynamic  control  of  tie-line  power  flow. 
The  redox  flow  batteries,  which  are  not  aged  to  the  frequent  charging  and  discharging,  have  a  quick 
response  and  outstanding  function  during  overload  conditions.  In  addition  to  leveling  load,  the  battery  is 
advantageous  for  secondary  control  in  the  power  system  and  maintenance  of  power  quality  of  distributes 
power  resources.  The  Bacterial  Foraging  Optimization  (BFO)  algorithm  is  use  to  optimize  the  parameters 
of  the  cost  functions  for  designing  the  integral  controller.  Compliance  with  North  American  Electric 
Reliability  Council  (NERC)  standards  for  Load  Frequency  Control  has  also  been  establishes  in  this  work. 
Simulation  studies  reveal  that  the  RFB  coordinate  with  IPFC  units  has  greater  potential  for  improving  the 
system  dynamic  performance. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  an  interconnected  power  system,  a  sudden  load  perturbation 
in  any  area  causes  the  deviation  of  frequencies  of  all  the  areas  and 
also  of  the  tie-line  powers.  This  has  to  be  corrected  to  ensure 
generation  and  distribution  of  electric  power  with  good  quality. 
This  is  achieved  by  Load  Frequency  Control,  also  known  as  Auto¬ 
matic  Generation  Control  (AGC).  The  main  objectives  of  AGC  [1,2] 
are  to  be  maintain  the  desired  megawatt  output  and  the 
nominal  frequency  in  an  interconnected  power  system  besides 
maintaining  the  net  interchange  of  power  between  control  areas  at 
predetermined  values.  The  paper  proposes  a  control  scheme  that 
ensures  reliability  and  quality  of  power  supply,  with  minimum 
transient  deviations  and  ensures  zero  steady  state  error.  The 
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importance  of  decentralize  controllers  for  multi  area  load- 
frequency  control  system,  where  in,  each  area  controller  uses 
only  the  local  states  for  feedback,  is  well  known.  The  stabilization 
of  frequency  oscillations  in  an  interconnected  power  system 
becomes  challenging  when  implements  in  the  future  competitive 
environment.  So  advanced  economic,  high  efficiency  and 
improved  control  schemes  [3,4]  are  required  to  ensure  the  power 
system  reliability.  The  conventional  load-frequency  controller  may 
no  longer  be  able  to  attenuate  the  large  frequency  oscillation  due 
to  the  slow  response  of  the  governor  [5].  The  recent  advances  in 
power  electronics  have  led  to  the  development  of  the  Flexible 
Alternating  Current  Transmission  Systems  (FACTS).  These  FACTS 
devices  are  capable  of  controlling  the  network  condition  in  a  very 
fast  manner  [6]  and  because  of  this  reason  the  usage  of  FACTS 
devices  are  more  apt  to  improve  the  stability  of  power  system. 
Several  FACTS  devices,  such  as  Thyristor  Controlled  Series  Capac¬ 
itor  (TCSC),  STATic  synchronous  COMpensator  (STATCOM), 
Thyristor  Controlled  Phase  Shifter  (TCPS),  Static  Synchronous 
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Series  Compensator  (SSSC),  Unified  Power  Flow  Controller  (UPFC), 
Interline  Power  Flow  Controller  (IPFC),  have  been  developed  in 
recent  decades  [7].  An  Interline  Power  Flow  Controller  consists  of 
a  set  of  converters  that  are  connected  in  series  with  different 
transmissions  lines  which  can  effectively  manage  the  power  flow 
control  in  multi  line  systems.  The  schematic  diagram  of  IPFC  is 
illustrated  in  Fig.  1.  The  Unified  Power  Flow  Controller  and  Inter¬ 
line  Power  Flow  Controller  consists  at  least  two  converters.  It  is 
found  that,  in  the  past,  much  effort  has  been  made  in  the  modeling 
of  the  UPFC  to  compensate  a  single  transmission  line  for  power 
flow  analysis  [8-10],  whereas  the  IPFC  is  conceived  for  the 
compensation  and  power  flow  management  of  multi-line  trans¬ 
mission  system.  Therefore,  UPFC  is  not  attractive  for  compensating 
multi-line  systems  from  economical  point  of  view  [11].  Interline 
Power  Flow  Controller  not  only  can  compensate  each  transmission 
line  separately  but  also  can  compensate  all  of  them  at  the  same 
time.  It  employs  a  number  of  Voltage  Source  Converters  (VSC) 
linked  at  the  same  DC  terminal;  each  of  them  can  provide  series 
compensation  for  its  own  line  [12].  In  this  way,  the  power  opti¬ 
mization  of  the  overall  system  can  be  obtain  in  the  form  of 
appropriate  power  transfer  through  the  common  DC  link  from 
over-loaded  lines  to  under-loaded  lines  [13].  A  simple  model  of 
IPFC  with  optimal  power  flow  control  method  to  solve  overload 
problem  and  the  power  flow  balance  for  the  minimum  cost  [14] 
has  been  proposed. 

In  this  paper  Interline  Power  Flow  Controller  is  being  install  in 
the  tie-line  between  any  interconnected  areas,  which  is  use  to 
stabilize  the  area  frequency  oscillations  by  high  speed  control  of 
tie-line  power  through  the  interconnections.  In  addition  it  can  also 
be  expect  that  the  high  speed  control  of  IPFC  can  be  coordinate  with 
slow  speed  control  of  governor  system  for  enhancing  stabilization 
of  area  frequency  oscillations  effectively.  Under  these  situations, 
the  governor  system  may  no  longer  be  able  to  absorb  the  frequency 
fluctuations.  In  order  to  compensate  for  sudden  load  changes,  an 
active  power  source  with  fast  response  such  as  redox  flow  batteries 
is  expect  as  the  most  effective  counter  measure.  The  redox  flow 
batteries  will,  in  addition  to  load  leveling,  a  function  conventionally 
assigned  to  them,  have  a  wide  range  of  applications  such  as  power 
quality  maintenance  for  decentralized  power  supplies.  The  redox 
flow  batteries  are  the  excellent  short-time  overload  output  and  the 
response  characteristics  possessed  in  the  particular  [15,16].  The 
effect  of  generation  control  and  the  absorption  of  power  fluctuation 
required  for  power  quality  maintenance  are  expected,  however  it 
will  be  difficult  to  locate  the  placement  of  RFB  alone  in  every 
possible  area  in  the  interconnected  system  due  to  the  economical 
reasons.  Therefore  RFB  coordinate  with  IPFC  are  capable  of 
controlling  the  network  conditions  in  a  very  fast  and  economical 
manner. 

Nowadays  power  system  complex  are  being  solved  with  the  use 
of  Evolutionary  Computation  (EC)  such  as  Differential  Evolution 
(DE)  [17],  Genetic  Algorithms  [GA],  Practical  Swarm  Optimizations 
[PSO]  [18],  Ant  Colony  Optimization  [ACO]  [19],  which  are  some  of 
the  heuristic  techniques  having  immense  capability  of  determining 


global  optimum.  Classical  approach  based  optimization  for 
controller  gains  is  a  trial  and  error  method  and  extremely  time 
consuming  when  several  parameters  have  to  be  optimized  simul¬ 
taneously  and  provides  suboptimal  result  [18,19].  Some  authors 
have  applied  GA  to  optimize  controller  gains  more  effectively  and 
efficiently  than  the  classical  approach.  But  the  premature  conver¬ 
gence  of  GA  degrades  its  search  capability.  Recent  research  has 
brought  out  some  deficiencies  in  using  GA,  PSO  based  techniques 
[20,21].  The  Bacterial  Foraging  Optimization  [BFO]  mimics  how 
bacteria  forage  over  a  landscape  of  nutrients  to  perform  parallel 
non  gradient  optimization  [22].  The  BFO  algorithm  is  a  computa¬ 
tional  intelligence  based  technique  that  is  not  large  affected  by  the 
size  and  non-linearity  of  the  problem  and  can  be  convergence  to 
the  optimal  solution  in  many  problems  where  most  analytical 
methods  faith  convergence.  A  more  recent  and  powerful  evolu¬ 
tionary  computational  technique  BFO  [23]  is  found  to  be  user 
friendly  and  is  adopted  for  simultaneous  optimization  of  several 
parameters  for  both  primary  and  secondary  control  loops  of  the 
governor.  In  this  study,  a  BFO  algorithm  is  used  to  optimizing  the 
integral  controller  gains  for  load  frequency  control  of  a  two  area 
thermal  power  system  without  and  with  IPFC  and  RFB.  To  obtain 
the  best  convergence  performance,  new  cost  function  is  derive  by 
using  the  tie-line  power  and  frequency  deviations  of  the  control 
areas  and  their  rates  of  changes  according  to  time  integral.  The 
main  function  of  LFC  is  to  regulate  a  signal  called  Area  Control  Error 
(ACE),  which  accounts  for  error  in  the  frequency  as  well  as  the 
errors  in  the  interchange  power  with  neighboring  areas.  Conven¬ 
tional  Load  Frequency  Control  uses  a  feedback  signal  that  is  either 
based  on  the  Integral  of  ACE  or  is  based  on  ACE  and  it’s  Integral. 
These  feedback  signals  are  used  to  maneuver  the  turbine  governor 
set  points  of  the  generators  so  that  the  generated  power  follows  the 
load  fluctuations,  however  continuously  tracking  load  fluctuations 
definitely  causes  wear  and  tear  on  governor’s  equipment,  shortens 
their  lifetime,  and  thus  requires  replacing  them,  which  can  be  very 
costly.  In  Ref.  [24]  deals  with  discrete-time  Automatic  Generation 
Control  of  an  interconnected  reheat  thermal  system  considering 
a  Area  Control  Error  New  (ACEN)  based  on  tie-power  deviation, 
frequency  deviation,  time  error  and  inadvertent  interchange.  This 
controller  can  effectively  regulate  time  error  (?)  and  inadvertent 
interchange  accumulations  (I).  Control  Performance  Criteria  (CPC) 
has  been  formerly  used  to  evaluate  AGC  performance.  The  Control 
Performance  Standard  (CPS)  is  specifically  designed  to  comply  with 
the  performance  standards  imposed  by  the  North  American  Electric 
Reliability  Council  for  equitable  operation  of  an  interconnected 
system.  Control  Performance  Standard  (CPS1)  and  Control  Perfor¬ 
mance  Standard  2  (CPS2)  are  derived  from  rigorous  theoretical 
basis.  CPS1  is  a  measurement  to  asses  the  performance  of  frequency 
control  in  each  area.  CPS2  is  designed  to  restrain  the  ACE  10-min 
average  value  and  in  doing  so  provides  a  means  to  limit  excessive 
unscheduled  power  flows  that  could  results  from  large  ACE.  In  this 
paper  a  novel  load  frequency  controller  is  present.  It  is  manipulate 
by  a  Fuzzy  logic  system  whose  rules  are  design  to  reduce  wear  and 
tear  of  the  equipment  and  assure  its  control  performance  is  in 
compliance  with  NERC’s  control  performance  standards  [25,26], 
CPS1  and  CPS2.  Considering  the  power  system  load  frequency 
control,  this  paper  establishes  a  fuzzy  logic  controller  to  predict  the 
future  frequency  of  the  target  object,  thus  forecasting  the  opti¬ 
mized  controller  is  designed,  which  follows  the  CPS  performance 
standards  through  the  fuzzy  logic  rules.  This  control  structure  is 
a  decentralized,  integral  type  controller  whose  parameter  is  auto¬ 
matically  tuned  using  Bacterial  Foraging  Optimization  algorithm. 
The  control  parameter  is  reduced  to  diminish  high  frequency 
movement  of  the  speed  governor’s  equipment  when  the  control 
area  has  high  compliance  with  NERC’s  standards.  When  the 
compliance  is  low,  the  control  parameter  is  raised  to  the  normal 
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value.  This  paper  adopts  CPS1  and  CPS2  as  input  to  the  fuzzy  logic 
controller  and  output  gain  of  fuzzy  controller  adjust  control 
parameter  gain  (Ki)  depending  up  on  NERC’s  compliance.  The 
simulation  results  show  that  the  dynamic  performance  of  the 
system  with  IPFC  and  RFB  is  improve  by  using  the  propose 
controller. 


the  converters  are  connected  together  via  a  common  dc  link  as 
shown  in  Fig.  1.  With  this  Inter-line  Power  Flow  Controller,  in 
addition  to  providing  series  reactive  compensation,  any  converter 
can  be  controlled  to  supply  real  power  to  the  common  dc  link  from 
its  own  transmission  line  [28].  It  can  effectively  manage  the  power 
flow  via  multi-line  transmission  system  by  the  injection  of  proper 
series  voltages  in  transmission  lines  with  the  aid  of  its  inverters. 


2.  Problem  formulation 

The  state  variable  equation  of  the  minimum  realization  model  of 
‘AT  area  interconnected  power  system  [27]  may  be  expressed  as 

x  =  Ax  +  Bu  +  rd 
y  =  cx 

where  x  =  [x\,  Ap^.x^^,  Ape(N_1)...xJ]T,  n  -  state  vector; 
n  =  f>  +  (N-l) 

Z  =  1 

u  =  [ui,  ...un]t  =  [APC1  ...Pcn]t,  N  -  control  input  vector; 

d  =  [d!,...dN]r=  [APD1  ...Pdn]t, 

N  -  disturbance  input  vector; 
y  =  [yi...yjv]T,  2N  -  measurable  output  vector; 

where  A  is  system  matrix,  B  is  the  input  distribution  matrix,  P  is  the 
disturbance  distribution  matrix,  C  is  the  control  output  distribution 
matrix,  x  is  the  state  vector,  u  is  the  control  vector  and  d  is  the 
disturbance  vector  consisting  of  load  changes.  In  order  to  ensure 
zero  steady  state  error  condition  an  integral  controller  may  suit¬ 
ability  designed  for  the  augmented  system.  To  incorporate  the 
integral  function  in  the  controller,  the  system  Eq.  (1)  is  augmented 
with  new  state  variables  defined  as  the  integral  of 
ACEj(/v|dt),  i  =  1,2, ...N.  The  augmented  system  of  the  order 
(N  +  n)  may  be  described  as 

x  =  Ax  +  Bu  +  Td  (2) 

where 


3.1.  Power  injection  model  of  inter-line  power  flow  controller 

In  this  section,  a  mathematical  model  for  IPFC  is  obtained  with 
the  power  injection  model  [29]  which  is  helpful  in  understanding 
the  impact  of  the  IPFC  in  the  power  system  during  the  steady  state. 
Furthermore,  the  IPFC  model  can  easily  be  incorporated  in  the 
power  flow  model.  Usually,  in  the  steady  state  analysis  of  power 
systems,  the  Voltage  Source  Converter  may  be  represented  as 
a  synchronous  voltage  source  injecting  an  almost  sinusoidal  voltage 
with  controllable  magnitude  and  angle.  Based  on  this,  the  equiva¬ 
lent  circuit  of  IPFC  is  developed  and  shown  in  Fig.  2.  Where  V/,  V, 
and  V  are  the  complex  bus  voltages  at  the  buses  i,  j  and  k  respec¬ 
tively,  Vsein  is  the  complex  controllable  series  injected  voltage 
source,  and  Zsein  (n  =  j,  k)  is  the  series  coupling  transformer 
impedance.  The  active  and  reactive  power  injections  at  each  bus 
can  be  easily  calculated  by  representing  IPFC  as  current  source.  For 
the  sake  of  simplicity,  the  resistance  of  the  transmission  lines  and 
the  series  coupling  transformers  are  neglected.  The  power  injec¬ 


tions  at  buses  are  summarized  as 

^mj,i  —  ^i^sein^insin(^i  —  ^sein)  (5) 

n  =jjc 

Qinj,/  =  —  VjVseinbincos(#/  —  #sein)  (6) 

n  =jjc 

Pinj,n  =  -VnVseinbmsm(6n  -  0sein)  (7) 

Qinj.n  =  VnVseinbincos(0n  -  0sein)  (8) 


X 


f  vdt 

)N  a- 

0 

c' 

B  = 

O' 

J 

L  *  J 

}n 

0 

A 

B 

and 


where 

Vx  =  VX/_6X  (x  =  i,j  and /<) 


The  problem  now  is  to  design  the  decentralized  feedback  control 
law 

Ui  =  -kjyi  i  =  1,2...JV  (3) 

The  control  law  in  Eq.  (3)  may  be  written  in-terms  of  v2-  as 

ui  =  -kiJvidt  i=l,2...N  (4) 

where  /q  is  the  integral  feedback  gain  vector,  v2  is  the  scalar  control 
output  of  area  i. 

3.  Operating  principle  of  inter-line  power  flow  controller 

In  its  general  form  the  Inter-line  Power  Flow  Controller  employs 
a  number  of  dc-to-ac  converters  each  providing  series  compensa¬ 
tion  for  a  different  line.  In  other  words,  the  Inter-line  Power  Flow 
Controller  comprises  a  number  of  Static  Synchronous  Series 
Compensators.  The  simplest  IPFC  consists  of  two  back-to-back  dc- 
to-ac  converters,  which  are  connects  in  series  with  two  transmission 
lines  through  series  coupling  transformers  and  the  dc  terminals  of 


Vsein  =  Vseinz0sein  (n  =  j,k). 

The  equivalent  power  injection  model  of  an  IPFC  is  shown  in 
Fig.  3.  As  Inter-line  Power  Flow  Controller  neither  absorbs  nor 
injects  active  power  with  respect  to  the  ac  system;  the  active  power 
exchange  between  the  converters  via  the  dc  link  is  zero,  i.e. 
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f/j 

K 

◄ — 

^tnlJ  +  JQ-WJ 

< - 

zseu 

K 

ftyj  +  jQb Vi 

+  jQtW* 

Fig.  3.  Power  injection  model  of  IPFC. 


Re(vsey/*+Vseifc4)  =  0  (9) 

where  the  superscript  *  denotes  the  complex  conjugate.  If  the 
resistances  of  series  transformers  are  neglected,  Eq.  (9)  can  be 
written  as 

E  Pinj,m  =  0  (10) 

m  =  ij,k 

Normally  in  the  steady  state  operation,  the  IPFC  is  used  to 
control  the  active  and  reactive  power  flows  in  the  transmission 
lines  in  which  it  is  placed.  The  active  and  reactive  power  flow 
control  constraints  are 


pm  -  p,r  =  o 

(ii) 

Qjii  -  Q,Tec  =  0 

(12) 

where  n  =  j,  k ;  P^ec,  Q^ec  are  the  specified  active  and  reactive 
power  flow  control  references  respectively,  and 

Pni  =  Re(vn4) 

(13) 

Qm  =  Im  (V„4) 

Thus,  the  power  balance  equations  are  as  follows 

(14) 

Pgm  +  Pinj,m  ~  Plm  +  ^line,??!  *®  0 

(15) 

Qgm  +  Qinj ,m  ~  Qlm  +  Qline,m  ®  6 

(16) 

where  Pgm  and  Qgm  <ire  generation  active  and  reactive  powers,  P/m 
and  Qim  are  load  active  and  reactive  powers.  Pime.m,  Qiine.m  are 
conventional  transmitted  active  and  reactive  powers  at  the  bus 
m  -  i,  j  and  k. 

4.  Operating  principle  of  redox  flow  batteries 

The  configuration  of  redox  flow  battery  is  shown  in  Fig.  4.  A 
sulphuric  acid  solution  containing  vanadium  ions  is  used  as  the 
positive  and  negative  electrolytes,  which  are  stored  in  respective 
tanks  and  circulated  to  battery  cell.  The  reactions  that  occur  in  the 
battery  cell  during  charging  and  discharging  can  be  express  simply 
with  the  following  equations  [15]. 


Fig.  4.  Principles  of  a  redox  flow  battery. 


.  charge  r 

Positive  Electrode  :  V4+  V5+  +  e-  (17) 

discharge 

charge  „ 

Negative  Electrode  :  V3+  +  e~  V2+  (18) 

discharge 

The  redox  flow  batteries  offer  the  following  features,  and  are 
suitable  for  high  capacity  systems  that  differ  from  conventional 
power  storage  batteries.  The  battery  reaction  only  involves 
a  change  in  the  valence  of  a  vanadium  ion  in  the  electrolyte.  There 
are  none  of  the  factors  which  reduce  the  battery  service  life  seen  in 
other  batteries  that  use  a  solid  active  substance,  such  as  loss  are 
electro  depositions  of  the  active  substance.  Further  more,  opera¬ 
tion  at  normal  temperatures  ensure  less  deterioration  of  the 
battery  materials  due  to  temperature.  The  system  configurations 
are  such  that  battery  output  (cell  section)  and  battery  capacity 
(tank  section)  can  be  separated,  therefore  the  layout  of  the  sections 
can  be  altered  according  to  the  place  of  installation.  For  example, 
the  tank  can  be  placed  underground.  The  design  can  be  easily 
modified  according  to  the  required  output  and  capacity.  The 
charged  electrolyte  is  stored  in  separate  positive  and  negative 
tanks  when  the  battery  has  been  charged,  therefore  no  self¬ 
discharge  occurs  during  prolonged  stoppage  nor  is  auxiliary 
power  required  during  stoppage.  Furthermore,  start-up  after  pro¬ 
longed  stoppage  requires  only  starting  of  the  pump,  thus  making 
start-up  possible  in  only  a  few  minutes.  The  electrolyte  (i.e.,  the 
active  substance)  is  sent  to  the  each  battery  cell  from  the  same 
tank,  therefore  the  charging  state  of  each  battery  cell  is  the  same, 
eliminating  the  need  for  special  operation  such  as  uniform 
charging.  So  that,  maintenance  is  also  easy  because  the  electrolyte 
is  relatively  safe  and  the  operating  are  at  normal  temperature  and 
assures  superior  environmental  safety.  Waste  vanadium  from 
generating  stations  can  be  used  so  it  can  be  superior  recyclability. 
Furthermore,  the  vanadium  in  the  electrolyte  can  be  used  semi¬ 
permanently. 

The  RFB  systems  are  incorporated  in  the  power  system  to 
suppress  the  load  frequency  control  problem  and  to  ensure  an 
improved  power  quality.  In  particular,  these  are  essential  for 
reusable  energy  generation  units,  such  as  wind  power  and  photo¬ 
voltaic  generator  units,  which  need  measures  for  absorption  of 
changes  in  output  and  to  control  flicker  and  momentary  voltage 
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drop.  With  the  excellent  short-time  overload  output  and  response 
characteristics  possessed  by  RFB  in  particular  [16],  the  effects  of 
generation  control  and  of  the  absorption  of  power  fluctuation 
needed  for  power  quality  maintenance  are  expected.  The  set  value 
of  the  RFB  has  to  be  restored  at  the  earliest,  after  a  load  disturbance 
so  that  the  RFB  unit  is  ready  to  act  for  the  next  load  disturbance.  The 
redox  flow  batteries  are  capable  of  ensuring  a  very  fast  response 
and  therefore,  hunting  due  to  a  delay  in  response  does  not  occur. 
For  this  reason,  the  AFi  was  used  directly  as  the  command  value  for 
LFC  to  control  the  output  of  RFB. 

5.  System  modeling  for  control  design 

The  redox  flow  battery  and  Inter-line  Power  Flow  Controller  are 
found  to  be  superior  to  the  governor  system  in  terms  of  the  faster 
response  against  the  frequency  fluctuations.  They  are  charged  with 
suppressing  the  peak  value  of  frequency  deviations  quickly  against 
the  sudden  load  change,  subsequently  the  input  to  the  governor 
system  are  updated  with  the  required  input  for  the  compensation 
of  the  steady  state  error  of  the  frequency  deviations.  Fig.  5  shows 
the  model  for  the  control  design  of  RFB  and  IPFC.  Where  the 
dynamics  of  governor  systems  are  eliminated  by  setting  the 
mechanical  inputs  to  be  constant  since  the  response  of  governor  is 
much  slower  than  that  of  RFB  or  IPFC.  The  redox  flow  battery  is 
modeled  as  an  active  power  source  to  areal  with  a  time  constant 
Trfb.  The  IPFC  is  modeled  as  a  tie  line  power  flow  controller  with 


a  time  constant  Tipfg  The  tie-line  power  modulated  by  the  IPFC 
flows  into  both  areas  simultaneously  with  different  signs  (+  and  -) 
since  the  responses  of  power  control  by  the  RFB  and  by  the  IPFC  are 
sufficiently  fast  compared  to  the  dynamics  of  the  frequency  devi¬ 
ations,  the  time  constants  Trfb  and  Tipfc  are  regarded  as  0  s  [5]  for 
the  control  design.  Then  the  state  equation  of  the  system  repre¬ 
sented  by  Fig.  5  becomes. 


'  4h  ' 

l 

O 

i 

HtS 

i 

1 

af2  J 

27 tT12  0  -27 tT12 

n  ankp2  1 

Tp 2  Tp2 

1 

fN 

kpi 

Tpi 

0 

0 


0 

aukp2 

h2 


APrfb 

APipfc 


(19) 


Here,  from  the  physical  view  point  it  is  noted  that  the  IPFC 
located  between  two  areas  is  effective  to  stabilize  the  inter-area 
oscillation  mode  only,  and  then  the  RFB  which  is  capable  of 
supplying  the  energy  into  the  power  system  should  be  suitable  for 
the  control  of  the  inertia  mode. 


Load  change  1 


Load  change  2 


Fig.  5.  Linearized  reduction  model  for  the  control  design. 
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5.1.  Control  design  of  redox  flow  batteries 


The  design  process  starts  from  the  reduction  of  two  area  system 
into  one  area  which  represents  the  Inertia  centre  mode  of  the  overall 
system.  The  controller  of  RFB  is  designed  in  the  equivalent  one  area 
system  to  reduce  the  frequency  deviation  of  inertia  centre.  The 
equivalent  system  is  derived  by  assuming  the  synchronizing  coef¬ 
ficient  T\2  to  be  large.  From  the  state  equation  of  APT12  in  Eq.  (19) 


APn2 

27rTi2 


Afr  -  af2 


(20) 


Setting  the  value  of  Ti2  in  Eq.  (20)  to  be  infinity  yields  AFi  =  AF2. 
Next,  by  multiplying  state  equation  of  A/7!  and  A F2  in  Eq.  (19)  by 
Tpi//<pi  and  Tp2/ai2  kP2  respectively,  then 


y^-A/7!  m  jT  AF"j  -  APn2  -  APipfc  +  APrfb  (21) 

Kpi  '<pi 

a^f^2  =  k^c bAF2  +  APt12  +  APlPFC  (22) 


Eq.  (19)  is  decomposed  into  two  decoupled  subsystems.  Where  the 
state  variable  APT12  is  duplicated  included  in  both  subsystems,  which 
is  the  reason  that  this  process  is  called  overlapping  decompositions. 
Then,  one  subsystem  which  preserves  the  inter-area  mode  is  repre¬ 
sented  by. 


'  AF,  ' 
APn2. 


1 


Tpi 

2ttF12 


-ftp  i 

hi 

0 


AF! 

APpiz 


+ 


-ftPi 

0 


[APipfc] 


(27) 


It  has  been  proved  that  the  stability  of  original  system  is  guar¬ 
anteed  by  stabilizing  every  subsystem.  Therefore  the  control 
scheme  of  IPFC  is  designed  to  enhance  the  stability  of  the  system 
Eq.  (27)  by  eigenvalue  assignment  method.  Flere  let  the  conjugate 
eigenvalue  pair  of  the  system  Eq.  (27)  be  a  ±j(3,  which  corresponds 
to  the  inter-area  mode.  The  control  purpose  of  the  IPFC  is  to  damp 
the  peak  value  of  frequency  deviation  in  area  1  after  a  sudden 
change  in  the  load  demand.  Since  the  system  Eq.  (27)  is  the  second 
order  oscillation  system,  the  percentage  overshoot  Mp  (new)  can  be 
specified  for  the  control  design.  Mp  (new)  is  given  as  a  function  of  the 
damping  ratio  by 


By  summing  Eqs.  (21)  and  (22)  and  using  the  above  relation 
AFi  =  AF2  =  AF 


The  load  change  in  this  system  APD  is  additionally  considered, 
where  C  is  constant,  here  the  control  APRFb  =  -I<rfb  AF  is  applied 
then. 


AF 


C 

S  +  A  +  Frfb^ 


APd 


(24) 


Where  A=  — ) 

\  fcpl  fcp2<3l2 ) 

B  =  T - - - T 

TP1  ,  TP2 

ftpl  ftp2a12 


TP1  ,  TP2  ^ 

kp  1  /<p2^12  J 


Since  the  control  purpose  of  RFB  is  to  suppress  the  deviation  of 
AF  quickly  against  the  sudden  change  of  APd,  the  percent  reduction 
of  the  final  value  after  applying  a  step  change  APD  can  be  given  as 
a  control  specification.  In  Eq.  (24)  the  final  values  with  JCrfb  =  0  and 
with  JCrfb  ^  0  are  C/A  and  C/(A  +  JCrfb#)  respectively  therefore  the 
percent  reduction  is  represented  by 


C/(A  +  /CrFB£)/(C/A)  =  ft/100 


(25) 


Mp(new)  =  e(  )  (28) 

The  real  and  imaginary  parts  of  eigenvalue  after  the  control  are 
expressed  by 


as  =  dwn  (29) 

0s  =  wn\f\  -  <52  (30) 

where  wn  is  the  undamped  natural  frequency,  by  specifying  Mp  and 
assuming  (3S  =  0,  the  desired  pair  of  eigenvalue  is  fixed.  As  a  result, 
the  eigenvalue  assignment  method  derives  to  feed  back  scheme  as 

APipfc  =  — /<i  AF1  — /<2APn2  (31) 

The  characteristic  polynomial  of  the  system  Eq.  (27)  with  state 
feedback,  which  is  given  by 

\XI  -  (A-  BI<)  |  =  0  (32) 

where  state  feedback  gain  matrix  I<  =  [/<i,/<2]-  The  desired  charac¬ 
teristic  polynomial  from  the  specified  eigenvalue  (^1,^2)  is  given  by 

(7-Mi)(A-m2)  =  0  (33) 

By  equating  the  coefficients  of  Eqs.  (32)  and  (33)  the  elements  k\,  /<2 
of  state  feedback  gain  matrix  I<  are  obtained. 

6.  Bacterial  foraging  optimization  technique 


For  a  given  R,  the  control  gain  of  RFB  is  calculated  as 


6.1.  Review  of  bacterial  foraging  optimization 


Krfb  =  4(10°-K)  (26) 


5.2.  Control  design  of  inter-line  power  flow  controller 

The  controller  for  the  IPFC  is  design  to  enhance  the  damping  of  the 
inter-area  mode.  In  order  to  extract  the  inter-area  mode  from  the 
system  Eq.  (19),  the  concept  of  overlapping  decompositions  is  applied. 
First,  the  state  variables  of  the  system  Eq.  (19)  are  classified  into  three 
groups,  i.e.  Xi  =  [AF^,  x2  =  [APn2],  x3  =  [AF2]  next,  the  system 


TheBFO  method  was  introduced  by  Possino  [22]  motivated  by 
the  natural  selection  which  tends  to  eliminate  the  animals  with 
poor  foraging  strategies  and  favor  those  having  successful  foraging 
strategies.  The  foraging  strategy  is  governed  by  four  processes 
namely  Chemotaxis,  Swarming,  Reproduction  and  Elimination  and 
Dispersal.  Chemotaxis  process  is  the  characteristics  of  movement  of 
bacteria  in  search  of  food  and  consists  of  two  processes  namely 
swimming  and  tumbling.  A  bacterium  is  said  to  be  swimming  if  it 
moves  in  a  predefined  direction,  and  tumbling  if  it  starts  moving  in 
an  altogether  different  direction.  To  represent  a  tumble,  a  unit 
length  random  direction  <p(j)  is  generated.  Let,  “j”  is  the  index  of 
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chemotactic  step,  “k"  is  reproduction  step  and  T  is  the  elimination 
dispersal  event.  6i(j,k,l ),  is  the  position  of  ith  bacteria  at  jth 
chemotactic  step  kth  reproduction  step  and  Ith  elimination 
dispersal  event.  The  position  of  the  bacteria  in  the  next  chemotatic 
step  after  a  tumble  is  given  by 

0f(/4-l,k,l)  =  0(j,k,l)  +  C(i)<p(j)  (34) 

If  the  health  of  the  bacteria  improves  after  the  tumble,  the 
bacteria  will  continue  to  swim  to  the  same  direction  for  the  spec¬ 
ified  steps  or  until  the  health  degrades.  Bacteria  exhibits  swarm 
behavior  i.e.  healthy  bacteria  try  to  attract  other  bacterium  so  that 
together  they  reach  the  desired  location  (solution  point)  more 
rapidly.  The  effect  of  swarming  [22]  is  to  make  the  bacteria 
congregate  into  groups  and  moves  as  concentric  patterns  with  high 
bacterial  density.  Mathematically  swarming  behavior  can  be 
modeled 


4.  Nc  —  the  number  of  iteration  to  be  undertaken  in  a  chemotactic 
loop  (Nc  >  Ns) 

5.  Nre  -  the  maximum  number  of  reproduction  to  be  undertaken. 

6.  Ned  —  the  maximum  number  of  elimination  and  dispersal 
events  to  be  imposed  over  bacteria 

7.  Ped  —  the  probability  with  which  the  elimination  and  dispersal 
will  continue. 

8.  The  location  of  each  bacterium  P  (1  -  p,  1  -  s,  1)  which  is 
specified  by  random  numbers  within  [-1,1] 

9.  The  value  of  C  (if  which  is  assumed  to  be  constant  in  our  case 
for  all  bacteria  to  simplify  the  design  strategy. 

10.  The  value  of  dattract»  ^attract*  hrepelent  Wrepelent*  It  is  to  be 
noted  here  that  the  value  of  dattract  and  /irepe lent  must  be  same 
so  that  the  penalty  imposed  on  the  cost  function  through  “Jcc” 
of  Eq.  (35)  well  be  “0”  when  all  the  bacteria  will  have  same 
value,  i.e.,  they  have  converged. 


Jcc(6,P(j,k,l))  =  jyccM(j.k.t)) 

i= 1 
S 

=  £ 

i=i 
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+£ 


dattract  exP(—Wattract)  ~  O'rn) 

m  =  l  7 


i  =  l 


^repelentexP  (  wrepelent) 


177=1 


(35) 


where  Jcc  ~  relative  distance  of  each  bacterium  from  the  fittest 
bacterium;  5  -  number  of  bacteria;  p  -  number  of  parameters  to  be 
optimized;  dm  -  position  of  the  fittest  bacteria;  dattract,  ^attract,  dre. 
peient.  wrepeient  —  different  co-efficients  representing  the  swarming 
behavior  of  the  bacteria  which  are  to  be  chosen  properly. 

In  Reproduction  step,  population  members  who  have  sufficient 
nutrients  will  reproduce  and  the  least  healthy  bacteria  will  die.  The 
healthier  population  replaces  unhealthy  bacteria  which  get  elimi¬ 
nated  owing  to  their  poorer  foraging  abilities.  This  makes  the 
population  of  bacteria  constant  in  the  evolution  process.  In  this 
process  a  sudden  unforeseen  event  may  drastically  alter  the 
evolution  and  may  cause  the  elimination  and/or  dispersion  to 
a  new  environment.  Elimination  and  dispersal  helps  in  reducing 
the  behavior  of  stagnation  i.e.,  being  trapped  in  a  premature  solu¬ 
tion  point  or  local  optima. 


6.2.  Bacterial  foraging  algorithm 

In  case  of  BFO  technique  each  bacterium  is  assigned  with  a  set  of 
variable  to  be  optimized  and  are  assigned  with  random  values  [A] 
within  the  universe  of  discourse  defined  through  upper  and  lower 
limit  between  which  the  optimum  value  is  likely  to  fall.  In  the 
proposed  method  integral  gain  K h  (i  =  1,  2)  scheduling,  each 
bacterium  is  allowed  to  take  all  possible  values  within  the  range 
and  the  cost  objective  function  which  is  represented  by  Eq.  (45)  is 
minimized.  In  this  study,  the  BFO  algorithm  reported  in 
Refs.  [22,23]  is  found  to  have  better  convergence  characteristics  and 
is  implemented  as  follows. 

6.2.1.  Step  1  -  initialization 

1.  Number  of  parameter  (p)  to  be  optimized. 

2.  Number  of  bacterial  (S)  to  be  used  for  searching  the  total 
region. 

3.  Swimming  length  (Ns),  after  which  tumbling  of  bacteria  will  be 
undertaken  in  a  chemotacticloop 


After  initialization  of  all  the  above  variables,  keeping  one  vari¬ 
able  changing  and  others  fixed  the  value  of  proposed  in  Eq.  (45) 
is  obtained  by  running  the  simulation  of  system  using  the  param¬ 
eter  contained  in  each  bacterium.  The  corresponding  to  the 
minimum  cost,  the  magnitude  of  the  changing  variable  is  selected. 
Similar  procedure  is  carried  out  for  other  variables  keeping  the 
already  optimized  one  unchanged.  In  this  way  all  the  variables  of 
step  1-  initialization  are  obtain  and  are  presented  below. 

S  =  6,  Nc  =  10,  Ns  *  3,  Nre  =  15,  Ned  =  2,  Ped 
—  0.25,  dattract  =  0.01,  Wattract  —  0.04,  hrepelent 
=  0.01,  andwrepelent  =  10.  p  =  2. 

6.2.2.  Step  2  -  iterative  algorithms  for  optimization 

This  section  models  the  bacterial  population  chemotaxis 
Swarming,  reproduction,  elimination,  and  dispersal  (initially, 
j  =  k  =  l  =  0)  for  the  algorithm  updating  0l  automatically  results  in 
updating  of  ‘P. 

1.  Elimination  —  dispersal  loop:  l  =  /  +  1 

2.  Reproduction  loop:  /<  =  /<  +  1 

3.  Chemotaxis  loop:  j  =  j  +  1 

a) .  For  i=l,  2...S,  calculate  cost  for  each  bacterium  i  as  follows. 

•  Compute  value  of  cost  J(ij,k,l ) 

Let  JswiijXl)  =  KijM)+JcJiOl(jXl)J(jM))  [i.e.,  add  on 
the  cell  to  cell  attractant  effect  obtained  through  Eq.  (35) 
for  swarming  behavior  to  obtain  the  cost  value  obtained 
through  Eq.  (45)]. 

•  Let  Jiast  =  JswiijXl )  to  save  this  value  since  a  better  cost 
via  a  run  be  found. 

•  End  of  for  loop. 

b) .  for  i  =  1,  2...S  take  the  tumbling/swimming  decision. 

•  Tumble:  generate  a  random  vector  A(i)e3ftp  with  each 
element  Am(i)  m  =  1,2,  ...p,  a  random  number  on  [-1,1]. 

•  Move  the  position  the  bacteria  in  the  next  chemotatic 
step  after  a  tumble  by  Eq.  (34)Fixed  step  size  in  the 
direction  of  tumble  for  bacterium  T  is  considered 

•  Compute  J(ij  +  1,/c,/)  and  then  let 

Jsw  ( ij + 1 5  k,  0  =J(iJ + 1  j  k,  l) 

+jcc(di(j+rk:i),P(j+rk.i))  (36) 


•  Swim: 

(i)  Let  m  =  0;  (counter  for  swim  length) 

(ii)  While  m  <  Ns  (have  not  climbed  down  too  long) 

•  Let  m  =  m  +  1 
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Table  1 

LFC  optimization  rules  based  on  CPS. 


Condition 


The  state  of  AGC  units 


CPS1  >  100%  and  CPS2  >  90% 
CPS1  <  100%  and  CPS2  >  90% 

CPS1  >  100%  and  CPS2  <  90% 
CPS1  <  100%  and  CPS2  <  90% 


No  optimization  adjusting 
ACE*AF  >  0  Optimization  adjusting 

ACE*AF  <0  No  optimization  adjusting 

Optimization  adjusting 
Optimization  adjusting 


'  If  JswO'j  +  l.fc./)  <  Just  (if  doing  better),  let 

Jiast  =Jsw{iJ  +  1 MJ)  and  let 


d'(j  +  \,k,l)  =  ei(j,k,l)+C(i) 


A(i) 


V/AV)A(0 


(37) 


where  C(z)  denotes  step  size;  A (z)  Random  vector;  AT(z)  Trans¬ 
pose  of  vector  A(z).  By  using  Eq.  (37)  to  compute  the  new 
J(ij  +  1,1c,/)  Else  let  zrz  =  Ns.  This  the  end  of  while  statement 

c).  Go  to  next  bacterium  (z  +  1 )  is  selected  if  z  =£  5  (i.e.  go  to  b) 
to  process  the  next  bacterium 

4.  If  j  <  Nc,  go  to  step  3.  In  this  case,  chemotaxis  is  continued  since 
the  life  of  the  bacteria  in  not  over 

5.  Reproduction. 

a)  For  the  given  k  and  l  for  each  z  =  1,2 ...S,  let 
^health  =  minje{i-Nc]{/sw(iJ,  k, /)}  be  the  health  of  the 
bacterium  z  (a  measure  of  how  many  nutrients  it  got  over  its 
lifetime  and  how  successful  it  was  in  avoiding  noxious 
substance).  Sort  bacteria  in  the  order  of  ascending  cost  Jheaith 
(higher  cost  means  lower  health). 

b)  The  Sr  =  S/2  bacteria  with  highest  Jheaith  values  die  and  other 
Sr  bacteria  with  the  best  value  split  [and  the  copies  that  are 
placed  at  the  same  location  as  their  parent]. 

6.  If  k  <  Nre,  go  to  2;  in  this  case,  as  the  number  of  specified 
reproduction  steps  have  not  been  reached,  so  the  next  gener¬ 
ation  in  the  chemotactic  loop  is  to  be  started 

7.  Elimination  —  dispersal:  for  z  =  1,  2...  S  with  probability  Pec i, 
eliminates  and  disperses  each  bacterium  [this  keeps  the 
number  of  bacteria  in  the  population  constant]  to  a  random 
location  on  the  optimization  domain. 

7.  North  American  electric  reliability  council’s  control 
performance  standards 


Fig.  7.  Membership  function  for  the  controller  outputs  (ai). 


and  it  is  defined  as  follows:  over  a  sliding  certain  period,  the 
average  of  the  “clock-min  averages”  of  a  control  area’s  Area  Control 
Error  divided  by  “10  times  its  area  frequency  bias”  times  the  cor¬ 
responding  “clock-  min  averages  of  the  interconnection  frequency 
error”  shall  be  less  than  the  square  of  a  given  constant  repre¬ 
senting  a  target  frequency  bound.  This  is  expressed  [32,33] 


AVG, 


period 


af, 


<q 


(38) 


where  ACE/:  Clock-average  of  ACE,  A Fz:  the  clock-  min  average  of 
frequency  deviation  from  schedule,  ft:  frequency  bias  of  the  control 
area,  e\  :  targeted  frequency  bound  for  CPS1,  z:  control  area  I  and  ( •  )i 
is  the  clock  1-  min  average. 

£i  is  a  constant  derived  the  historical  frequency  record  of 
a  control  area.  It  is  the  root  mean  square  of  one-  min  average 
frequency  deviation  from  a  schedule  based  on  frequency  perfor¬ 
mance  over  an  averaging  period  of  a  year.  The  period  (rz)  is  defined 
as  one  year  for  control  area  evaluation  or  one  month  for  the  report 
of  NERC.  To  calculate  CPS1  (/Ccpsi),  a  dimensionless  compliance 
factor  (/<cf)  is  defined  as: 


KCf 


ACE, 


AFi 


1/  1 


rzcf 


CPS1  is  then  obtained  from  the  following  equation: 

/Ccpsi  =  (2-/Ccf)*100% 


(39) 


(40) 


In  1997,  the  North  American  Electric  Reliability  Council  (NERC) 
proposed  new  control  performance  standards  [25,26]  CPS1  and 
CPS2  to  evaluate  the  control  area  performance  in  normal  inter¬ 
connected  power  system  operation.  Each  control  area  is  required  to 
monitor  its  control  performance  and  report  its  compliance  CPS1 
and  CPS2  to  NERC  [30,31]  at  the  end  of  each  month. 

7.1.  Control  performance  standard  l 

Control  Performance  Standard  1  assesses  the  impact  of  Area 
Control  Error  on  frequency  over  a  certain  period  window  or  horizon 


The  fundamental  requirement  for  CPS1  is  that  performance,  as 
measured  by  percentage  compliance  must  be  at  least  100%. 

1.  When  /Ccpsi  >  200%,  which  means  JCCf  <  0,  there 

(ACEi*AFi)  <0.  Under  this  condition,  ACE  facilitates  the 
frequency  quality. 

2.  When  100%  <  /Ccpsi  <  200%,  which  means  0  <  /Ccf  <  1,  there  is 
0  <  J][(ACE//  -  10ft)1*AF1]  <  nf\.  The  Control  Performance 
Standard  1  standard  is  satisfied. 

3.  When  /Ccpsi  <  100%,  which  means  /Ccf  >  h  there  is 
X)[(ACEf/  -  lOft^  *AF!]  >  nf\.  ACE  has  exceeded  the  permitted 


Fig.  6.  Membership  function  for  the  input  variables  (CPS1,  CPS2). 
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Table  2 

Fuzzy  logic  rules. 


CPS2 

CPS1 

PS 

PM 

PB 

PS 

ZE 

PS 

PS 

PM 

PS 

PB 

PM 

PB 

PB 

PVB 

PVB 

range  so  that  it  has  a  bad  effect  on  the  frequency  and  quality  of 
power  grid. 

7.2.  Control  performance  Standard  2 

The  Control  Performance  Standard  2  (CPS2),  limits  the  magni¬ 
tude  of  short-term  ACE  values.  It  requires  the  10-min  averages  of 
a  control  area  s  ACE  be  less  than  a  given  constant  (Lio),  as  in  the 
equation  below: 

AVG10min(ACEl)  <  L10  (41) 

where,  L10  =  1.65r10\/(-10,ftj)(-10/y.  Note  that  j 8S  is  the 
summation  of  the  frequency  bias  settings  of  all  control  areas  in  the 
considered  interconnection,  and  £i0  is  the  target  frequency  bound 
for  CPS2.  To  comply  with  this  standard,  each  control  area  must  have 
its  compliance  no  less  than  90%.  A  compliance  percentage  is 
calculated  from  the  following  equation: 

AVGiomin(ACEj)  //p\ 

kCPS2  S3  - - -  [42 ) 

MO 

In  order  to  meet  the  requirements  of  the  power  grid  frequency 
quality,  the  average  ACE  value  during  10-min  in  each  control  region 
should  be  in  the  normal  distribution  as: 


Table  3 

RFB  and  IPFC  control  design  results. 


IPFC 

RFB 

1.  Eigenvalue  of  system  Eq.  (27) 

71  =  -0.25  +  jl.8081 

72  =  -0.25  -jl.8081 

2.  Inter-area  mode,  Mp  =  63.25% 

3.  Design  specification,  MP(new)  =  10% 

4.  New  Eigenvalue  of  the  system  Eq.  (27) 

1.  Design  specification, 

R  =  2.4 

2.  RFB  gain  by  Eq.  (26), 
Krfb  =  0.67 

71  (new)  =  -1.323  +  jl.8081 

72(new)  =  -1.323  -  jl.8081 

5.  State  Feed  back  Gain  Eq.  (31) 

[IU,I<2 ]  =  [-0.432,  -0.532] 

a  =  £lo^/(-10ft)(-10  ft) 

(43) 

7.3.  Optimization  rules  based  on  control  performance  standards 

Suppose  Control  Performance  Standard  1  >  100%  and  Control 
Performance  Standard  2  >  90%  to  be  goal  of  the  LFC  control 
strategy.  Table  1  shows  LFC  optimization  rules  based  on  Control 
Performance  Standard. 

7.4.  Fuzzy  logic  design 

Fuzzy  logic  rules  are  designed  to  manipulate  the  conventional 
integral-  type  load  frequency  control  to  achieve  two  objectives:  (i) 
minimize  equipments’s  wear  and  tear  and  (ii)  comply  with  NERC,  s 
Control  Performance  Standard  1  and  Control  Performance  Standard 
2.  The  control  structure  for  each  area  is 


I  * 

-Ki 

X 

J  ► 

w 

s 

controller  2 


1  * 

1 

1 

Kr3Ti3.s+1 

\7 

J  ► 

Tg3.s+1 

w 

Tt3.s+1 

W 

Tr3.s+1 

governor  3 

turbine  3 

reheat  3 

1 

1 

Kr4Ti4.s+1 

Tg4.s+1 

_ 

Tt4.s+1 

Tr4.s+1 

governor  4 


Kps  2 
Tps2.s+1 


powersystem  2 


Fig.  8.  Transfer  function  model  of  two-area  interconnected  thermal  reheat  power  system  with  IPFC  and  RFB  units. 
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Table  4 

Tuned  parameter  of  the  control  system  for  three  case  studies. 


Two  area  interconnected 
power  system 

Feedback 
gain  (I<i) 

Cost  function 
value  (J) 

Case:  1 

Integral  controller  only 

Case:  2 

1.0848 

0.5879 

Integral  controller  with 

IPFC  unit  only 

Case:  3 

1.1058 

0.4627 

Integral  controller  with 

IPFC  and  RFB  units 

1.6321 

0.0146 

Uj  =  APci  =  MjKji  J  ACEj  dt  (44) 

where  A Pc  is  the  governor  set  point  or  raise/lower  signal,  I<\  the 
integral-control  parameter  and  a  is  set  using  fuzzy  logic  and  called 
fuzzy  gain.  This  paper  uses  information  that  reflects  compliance 
with  CPS1  and  CPS2  as  the  input  to  the  fuzzy  rules.  The  proposed 
fuzzy  logic  will  lower  the  integral  gain  I<\  when  the  control  area  has 
high  compliance.  On  the  other  hand,  that  the  integral  gain  I<\  will  be 


x  io-3 


Fig.  9.  a.  Dynamic  responses  of  the  frequency  deviations  of  area  1.  b.  Dynamic 
responses  of  the  frequency  deviations  of  area  2.  c.  Dynamic  responses  of  the  tie-line 
power  deviations. 


Fig.  10.  a.  Dynamic  responses  of  the  control  input  deviations  of  area  1.  b.  Dynamic 
responses  of  the  control  input  deviations  of  area  2. 

increased  when  the  compliance  with  CPS1  of  the  control  area  is 
low.  According  to  the  optimized  rules  from  Table  1,  the  membership 
functions  of  CPS1,  CPS2  and  a\  could  be  defined  as  Figs.  6  and  7. 
Fuzzy  rules  are  summarized  in  Table  2. 


Time  (s) 

Fig.  11.  a.  Dynamic  responses  of  the  required  additional  mechanical  power  generation 
of  area  1.  b.  Dynamic  responses  of  the  required  additional  mechanical  power  gener¬ 
ation  of  area  2. 
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8.  Simulations  results  and  observations 

The  effect  of  the  propose  redox  flow  batteries  coordinate  with 
Interline  Power  Flow  Controller  on  Load  Frequency  Control  inves¬ 
tigate  in  the  two  equal  area  interconnected  thermal  power  system, 
In  each  area  consists  of  two  reheat  units  as  shown  in  Fig.  8.  In  this 
study  the  active  power  model  of  IPFC  controllers  is  fitted  in  the  tie¬ 
line  near  areal  to  examine  its  effect  on  the  power  system  perfor¬ 
mance.  Then  RFB  unit  is  installing  in  area  1  and  coordinate  with 
IPFC  controller  for  LFC  to  study  its  performance  of  system.  The 
nominal  parameters  are  given  in  Appendix  A.  The  control  param¬ 
eters  of  redox  flow  batteries  and  Interline  Power  Flow  Controller 
are  design  as  shown  in  Table  3.  The  following  object  of  obtaining 
the  optimal  solutions  of  control  inputs  is  taken  an  optimization 
problem,  and  the  novel  cost  function  [34]  in  Eq.  (45)  are  derive  by 
using  the  frequency  deviations  of  control  areas  and  tie  line  power 
changes.  The  integral  controller  gain  is  tune  with  BFO  algorithm  by 
optimizing  the  solutions  of  control  inputs.  The  simulations  are 
realized  for  1%  step  load  perturbation  in  areal.  The  results  are  ob¬ 
tained  by  MATLAB  7.01  software  and  50  iterations  are  chosen  for 
converging  to  solution  in  the  BFO  algorithm. 


Time  (s) 


Fig.  12.  a.  Dynamic  responses  of  the  frequency  deviations  of  area  1  considering  CPS1 
and  CPS2.  b.  Dynamic  responses  of  the  frequency  deviations  of  area  2  considering  CPS1 
and  CPS2.  c.  Dynamic  responses  of  the  tie  line  power  deviation  considering  CPS1  and 
CPS2. 


T 

J  =  |{(ftA/i)2+(foA/2)2+(APtiei2)2}  (45) 

0 

The  BFO  algorithm  optimal  gains  of  Integral  controllers  (Kji  Kn) 
for  three  case  studies  are  determining  using  BFO  algorithm.  From 
the  simulation,  the  tuned  parameters  of  the  control  system  are 
found  and  shown  in  Table  4.  These  controllers  are  implementing  in 
an  interconnected  two-area  power  system  without  IPFC  and  RFB, 
with  IPFC  only,  IPFC  coordinate  RFB.  Relative  compliance  of  the 
propose  controller  based  LFC  schemes  to  the  NERC  standards  haven 
been  establish  for  the  above  power  system.  In  the  present  work, 
variation  of  load  in  only  areal  has  been  considered.  The  compliance 
factor  (J<cf)  is  compute  for  1-h  in  Eq.  (39).  Assuming  that  the 
response  of  the  controller  to  the  load  variations  for  the  year  will  be 
similar  to  that  obtained  during  the  sample  period  of  1-h.  According 
to  the  compliance  factor  (J<cf)  value  to  calculate  Control  Perfor¬ 
mance  Standard  1,  as  defined  in  Eq.  (40).  In  addition,  the  Control 
Performance  Standard  2,  as  defined  in  Eq.  (42),  has  also  been 
computed.  CPS1  and  CPS2  have  two  input  of  fuzzy  logic  controller, 
fuzzy  logic  rules  are  design  according  to  their  Compliance  the 
optimum  integral  gain  Kj  changes  depends  up  on  fuzzy  gain  (a).  It  is 
seen  from  Figs.9-11;  it  is  evident  that  the  dynamic  responses  have 
improved  significantly  with  the  use  of  information  that  reflects 
compliance  with  CPS1  and  CPS2.  This  algorithm  will  significantly 
reduce  wear  and  tear  of  the  equipment  since  movements  of  the 
governor  set  point  or  raise/lower  signal  (APC)  generate  from  the 
integral  controller  are  less  frequent  when  the  control  area  has  high 
compliance  or  that  values  of  1-min  average  compliance  factor  (CF) 
or  accumulatively  average  compliance  factor  (CF  ac)  is  less  than 
unity.  Simulation  results  are  shown  in  Figs.  12-14;  it  can  be 
observed  that  the  oscillations  in  area  frequencies  and  tie-line 
power  deviations  have  decreases  to  a  considerable  extent  as 
compare  to  that  of  the  system  without  IPFC  controllers.  Moreover 
the  RFB  locate  in  area  1  which  as  coordinate  with  IPFC  controllers. 


Fig.  13.  a.  Dynamic  responses  of  the  control  input  deviations  of  areal  considering  CPS1 
and  CPS2.  b.  Dynamic  responses  of  the  control  input  deviations  of  area2  considering 
CPS1  and  CPS2. 
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Time  (s) 


Time  (s) 

Fig.  14.  a.  Dynamic  responses  of  the  required  additional  mechanical  power  generation 
of  areal  considering  CPS1  and  CPS2.  b.  Dynamic  responses  of  the  required  additional 
mechanical  power  generation  of  area  2  considering  CPS1  and  CPS2. 

The  gain  of  RFB  is  calculating  using  Eq.  (26)  for  given  value  of  speed 
regulation  coefficient  (ft)  as  shown  in  Table  3.  Redox  flow  batteries 
improving  the  inertia  mode  oscillation.  The  settling  times  and  peak 
over/under  shoot  for  the  frequency  deviations  in  each  area  and 
tie-line  power  deviations  for  three  case  studies  are  tabulated  in 


Table  5 

Comparison  of  the  system  performance  for  the  four  case  studies. 


Two  area  interconnected  Setting  time  (is 

)  in(s) 

Peak  over/under  shoot 

power  system 

A  Fi 

A  F2 

APtie 

A Pi  (Hz)  A F2  (Hz)  APtie 

(p.u.MW) 

Case:  1 

Integral 

16.37 

15.19 

16.48 

0.0226 

0.0185 

0.0064 

controller 

only 

Case:  2 

Integral 

13.76 

12.72 

12.44 

0.0210 

0.0171 

0.0061 

controller 

considering 

CPS1 

and  CPS2 

Case:  3 

Integral 

8.255 

7.882 

11.66 

0.0126 

0.0108 

0.0039 

controller 

considering 

CPS1  and 

CPS2  with  IPFC 
units  only 

Case:  4 

Integral 

5.506 

7.727 

7.625 

0.0059 

0.0055 

0.0022 

controller 

considering 

CPS1  and 

CPS2  with 

IPFC 

and  RFB  units 

Table  5.  Fig.  14  shows  the  generation  responses  considering 
apfn  =  apfi2  =  0.5  and  apf2i  =  apf22  =  0.5  for  three  case  studies  and 
RFB  having  A/j  as  the  control  logic  signals.  As  the  load  disturbances 
have  occurred  in  area  1,  at  steady  state,  the  powers  generated 
by  generating  units  in  both  areas  are  in  proportion  to  the 
area  participation  factors.  Therefore,  at  steady  state 
APgiss  m  APg2ss  =  APDi*apfn  =  0.01*0.5  =  0.005  pu.MW  and 
similarly,  APg3ss  =  APg4ss  =  APD2  *  apf2i=0  *  0.5=0  pu.MW.  From 
Tables  4  and  5,  it  can  be  observe  that  the  controller  design  using 
BFO  algorithm  for  two  area  thermal  reheat  power  system  with  RFB 
coordinate  with  IPFC  have  not  only  reduces  the  cost  function  but 
also  ensure  better  stability,  as  they  possesses  less  over/under  shoot 
and  faster  settling  time  when  compare  with  the  output  response  of 
the  system  used  controller  design  for  the  two  area-two  unit 
thermal  reheat  power  system  without  IPFC  and  RFB  units.  More¬ 
over  the  RFB  units,  suppresses  the  peak  frequency  deviations  of 
both  areas,  and  continue  to  eliminate  the  steady  state  error  of 
frequency  deviations. 


9.  Conclusions 

This  paper  presents  a  design  of  load  frequency  control  for 
interconnected  power  systems  using  Bacterial  Foraging  Optimiza¬ 
tion  algorithm.  This  algorithm  is  employ  to  achieve  the  optimal 
parameters  of  the  compensating  units.  The  design  objectives  are  (i) 
to  comply  with  the  North  American  Electric  Reliability  council’s 
control  performance  standards;  CPS1  and  CPS2,  (ii)  to  reduce  wear 
and  tear  of  generating  unit’s  equipments,  and  (iii)  to  design  feasible 
control  structure.  The  control  structure  is  select  to  be  consistent 
with  conventional  LFC  designs  and  it  is  chosen  to  be  a  feedback  of 
a  signal  proportional  to  the  integral  of  the  Area  Control  Error.  The 
gain  of  this  integral-type  controller  consists  of  products  of  two 
terms,  a  conventional  control  gain  and  fuzzy  gain.  The  fuzzy  gain  is 
set  using  fuzzy  logic  rules.  Which  are  developed  to  comply  with 
NERC’s  standards  and  to  manipulate  the  generator’s  set  points  only 
if  need  be  to  reduce  excessive  maneuvering  and  hence  minimize 
the  cost  of  operation  and  maintenance  associate  with  LFC.  A 
sophisticated  Load  Frequency  Control  by  redox  flow  batteries 
coordinate  with  Interline  Power  Flow  Controller  has  proposed  for 
a  two  area  interconnected  Reheat  Thermal  Power  System.  This  BFO 
algorithm  has  faster  converging  algorithm  and  will  reduce 
computational  burden.  A  control  strategy  has  been  proposed  to 
adjust  the  voltage  of  IPFC  which  in  turn  controls  the  inter-area  tie¬ 
line  power  flow.  Simulation  results  reveal  that  the  first  peak 
frequency  deviation  of  both  areas  and  tie-line  power  oscillations 
following  sudden  load  disturbances  in  either  of  the  areas  can  be 
suppress  a  controlling  the  series  voltage  of  IPFC.  Moreover,  the  tie¬ 
line  power  flow  control  by  an  IPFC  unit  has  to  be  efficient  and 
effective  for  improving  the  dynamic  performance  of  load  frequency 
control  of  inter  connected  power  system  than  that  of  the  system 
without  IPFC  controllers.  The  advantages  of  the  expect  RFB  over 
existing  power  system  in  the  LFC  applications  were  examined.  RFB 
contributes  a  lot  in  promoting  the  efficiency  of  overall  generation 
control  through  the  effect  of  the  use  in  load  leveling  and  the 
assurance  of  LFC  capacity  after  overload  characteristic  and  quick 
responsiveness.  Simulation  results  reveal  that  the  design  concept  of 
damping  the  inertia  mode  and  inter-area  mode,  the  co-ordinate 
control  is  effective  to  suppress  the  frequency  deviation  of  two 
area  system  simultaneously.  It  may  be  therefore  be  conclude  that, 
the  redox  flow  batteries  with  a  sufficient  margin  of  LFC  capacity 
absorbs  the  speed  governor  capability  in  excess  of  falling  short  of 
the  frequency  bias  value  and  tie-line  power  flow  control  by  an  Inter 
line  Power  Flow  Controller  be  expect  to  be  utilize  as  a  new  ancillary 
service  for  the  stabilization  of  the  tie-line  power  oscillations  even 
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under  the  congestion  management  environment  of  the  power 
system. 
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Appendix  A 

(A)  Data  for  thermal  reheat  power  system  [27,32,33]. 

Rating  of  each  area  =  2000  MW,  Base  power  =  2000  MVA, 
fo  =  60  Hz,  Ri  =  R2  =  R3  =  Ra  =  2.4  Hz/p.u.MW, 
Tgi  =  Ig2  =  Tg3  =  Tg4  =  0.08  s,  Tri  =  Tr 2  =  Tr \  —  Tr 2  10  s, 
Tti  =  Tt2  =  Tts  =  Tt4  =  0.3  s,  /Cpi  =  I<p2  =  120  Hz/p.u.MW, 
Tpi  =  rp2  =  20  s,  ft  =  fo  «  0-425  p.u.MW/Hz, 

JCri  =  /Cr2  =  /Cr3  =  /Cr4  =  0.5,  27tTi2=0.545  p.u.MW/Hz,  a12  =  -1, 
APDi  =  0.01  p.u.MW,  £1  =  18  mHz,  ci0  =  5.7  mHz. 

(B)  Data  for  IPFC  [12]:  Tipfc  =  0.01  s; 

(C)  Data  for  RFB  [16]:  TrFb  =  0  s. 
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